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We consider quantum devices for turning a finite number TV of d-level quantum systems in the 
same unknown pure state o into M > N systems of the same kind, in an approximation of the M- 
fold tensor product of the state a. In a previous paper it was shown that this problem has a unique 
optimal solution, when the quality of the output is judged by arbitrary measurements, involving 
also the correlations between the clones. We show in this paper, that if the quality judgement is 
based solely on measurements of single output clones, there is again a unique optimal cloning device, 
which coincides with the one found previously. 



According to the well known "no-cloning theorem" [pj perfect copying of quantum information is impossible, i.e. 
there is no machine which takes a quantum system as input and produces two systems of the same kind, both of them 
indistinguishable from the input. However, from the point of view of practical applications in Quantum Information 
Theory this Theorem by itself is not very useful, because it only asserts that the cloning task cannot be performed 
exactly — but then no task can be performed exactly by real devices. The fundamental importance of the No-Cloning 
Theorem is expressed much better in stronger versions of the Theorem, which also give explicit lower bounds on the 
error made in any attempt to build a cloning device. Some such bounds have been established, as well Even 
more insight into the cloning problem is given by results showing how to minimize the error, i.e., how to construct 
optimal cloning devices (l|-§| ■ Other recent related work can be found in [p|-|l5| . 

In this paper we consider cloning devices, which take as input a certain number N of identically prepared systems, 
and produce a larger number M of systems as output. Again, the cloning task is to make the output state resemble 
as much as possible a state of M systems all prepared in the same state as the inputs. This variant of the problem is 
of interest as a "quantum amplifier" . It also has a better chance of reasonable success than a cloning device operating 
on single input systems: in the limit of many input systems the device can make a good statistical estimate of the 
input density matrix and hence produce arbitrarily good clones. 

Different variants of this problem arise by different choices of the type of systems and the set of states which should 
be copied, e.g. pure vs. mixed states, or a finite number of states arising in a cryptographic protocol. In the present 
paper we are exclusively concerned with the cloning of arbitrary unknown pure states. 

A second choice to be made is the precise notion of approximation between the output states of the cloning 
device and the (inattainable) target state. Apart from technicalities the basic issue here is whether the full states 
are compared, or only the one-clone marginals. Approximation in the first sense means that the expectations of all 
observables, including those testing correlations and entanglement between different clones, are close in the two states 
being compared. On the other hand, approximation in the second sense means closeness of expectations of single clone 
observables only. Perhaps this second condition has more of the flavour of the No-Cloning Theorem, since in that 
theorem, too, the requirement is that each single (!) clone be indistinguishable from the input. 

In p6[ we showed that the pure state cloning problem with all-particle test criterion has a unique optimal solution. 
In this paper we show the same for the single-particle test criterion, and that the two optimal cloning devices are 
actually the same. The difference between the two results may not seem great. However, the result in the present 
paper required much heavier mathematical machinery, and we believe it to be considerably deeper. The reason is 
that one-particle tests by far do not exhaust the linear space of M -particle observables. In particular, all correlations 
of the doner's output are ignored by the test, which would make a unique optimal solution appear rather unlikely. 
Nevertheless, this is what we prove. 



I. INTRODUCTION 



f 



II. STATEMENT OF THE PROBLEM AND MAIN RESULT 



Let us start with a precise formulation of the question we are going to consider. First of all, we will study throughout 
this paper only d-level systems with arbitrary but finite d. Hence the one particle Hilbert space TL we are using is 
TL = C d . The Hilbert space for the input to the cloning device is therefore the TV- fold tensor product TL® N of TL with 
itself. In fact, because we only consider tensor powers of pure states as inputs, it suffices to take the subspace of TL® N 
spanned by vectors of the form tp® N with ip G TL. This is precisely the "Bose" subspace TL® N C TL® N , i.e., the space 
of vectors invariant under all permutations. The output Hilbert space will be TL® M with M > N. On this space we 
cannot impose an a priori symmetry restriction, although such a restriction will come out a posteriori, as a special 
property of optimal cloning devices. 

A cloning map is a completely positive, unital map T : B(TL® M ) -> B(Hf N ). This describes the action of the device 
on observables. Its (pre-)dual, describing the same operation in terms of states, will be denoted]^] by T* : B*(TL® N ) — » 
B*(H® M ). If we identify states with density operators, this means that tr(pT(^4)) = tr(T*(p)A) for arbitrary density 
operators p and observables A. The input of the cloning device are N systems, prepared independently according to 
the same state a. Thus the overall input state is a® N . We will assume a to be pure, i.e., the density matrix of a is 
a one-dimensional projection onto a wave vector ip S TL, say. Then a® N is the projection onto the vector ip® N . The 
output of the cloning device is the state T*(a® N ), which is a (generally entangled) state of M > N systems. Our aim 
is to design T so that the output states T*(<7® ) approximate the product states <j® M . 

The one particle observables, on which the comparison will be based, will be written as a iff) = ® a ® 

j<a(M-k) g b(H® m ), for all a e B{H). Thus the optimal cloning problem for arbitrary pure input states T is to make 
the expectations 

tv(a (k )T*(*® N )) = tr(T(a {k) )a® N ) = ,T(a {k) )^ N ) 
and tr(a(fc)cr®^ / ) = tr(acr) = (ip,aip) 

as similar as possible for arbitrary one-particle observables a and one-particle vectors ip. Of course, when taking a 
supremum over such differences, the size of a has to be constrained somehow. We will choose the constraint < a < 1, 
which is to say that the above two expressions have an immediate interpretation as probabilities. The largest difference 
of such probabilities is now the error functional for cloning maps, which we will seek to minimize: 

A onc (T) = sup \(ip® N ,T(a [k) )iP® N ) - (ip,aip)\ (1) 

where the supremum is taken over all ip £ 7i with ||^|| = 1, all operators a <E B(H) with < a < I, and all integers 
1 < k < M. 

The corresponding quantity based on tests of the full output state (including correlations) is 

A all (T) = sup sup |tr (T(A)a® N ) - ix(Aa 9M )\ , 

A fJ,pure 

where the supremum is taken over all A £ B(H® M ) with < A < 1 and over all pure states a E £>*(7i). Due to the 
properties of the trace norm || • ||i this functional can be expressed by 

A all (T)= sup \\T.{a 9N ) - <r® M ||!. (2) 

CT,purc 

It turns out that there is exactly one cloning map T which minimizes this error functional. This can be proven with 
a minor adaptation of the arguments in flfj], which start from a slightly different criterion, namely the maximization 
of the "fidelity" T{T*) = su P(T purc (cr® M T»(cr (8Ar )). The unique solution T = f minimizing (|), or maximizing .F(T*), 
is best expressed in terms of its action on states, i.e., 

T*(p) = ^ ] S M ( P ®l M - N )S M . (3) 



1 In [|16| the symbol T was used for T*. In contrast to the key arguments in the present paper are phrased more readily 
in terms of the map on observables than in terms of the map on states. Therefore, we decided to change this notation, which 
then also agrees better with the usage for completely positive maps on operator algebras. 
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Here d[N] = ( d+ ^ 1 ) denotes the dimension of the symmetric subspace TiS_ , Sm is the projection from Ti® M to 

7i® M , and p is an arbitrary density operator on Ti® N . In |l6| we also computed the one-site restriction of the output 
states of this doner: 



tr 



(T(a k )a® N 



7(T)a(a) + (l- 7 (T))tr(a)/d , 



where 



7cn 



N M + d 



N + d M 



is the so-called Black Cow factor of T, interpreted as a "shrinking factor of the Poincare sphere" in the discussions of 
the qubit (d — 2) case. This makes it easy to verify the case of equality in the following Theorem, which is our main 
result. 

1 Theorem. For any cloning map T : B(H® M ) -> B{Tif N ) we have 



A onc (T) > 



d- 1 



1 - 



N M + d 



N + d M 



with equality iffT — T with T from equation (^). 



III. FINDING THE OPTIMAL CLONING MAP 



A. Reduction to the covariant case 



In this section we will give the proof of our main theorem, apart from some group theoretical Lemmas, which will 
be proved in Appendix A. Throughout, the symmetry of sitewise unitary rotation of clones and input states will play 
a crucial role. The necessary background information on unitary representations of SU(d) will also be supplied in 
Appendix A. 

We establish some notation first. By U(d) we will denote the group of unitary d x d-matrices, i.e., the unitary 
group on our underlying one-particle space Ti = C . Unitary representations of this group will be denoted by the 
letter 7r with suitable indices. 7r n is the defining representation on C d , and its n th tensor power, acting on 7^®^ by 
the operators ir^ N (u) — u® N is 7rjf N . The restriction of this representation to the symmetric subspace 7i^. N will be 
denoted by tt+. Thus a cloning map T : B{H® M ) -> B{H% N ) is called U '(d)- covariant, if 

T(7r® M (u)An® M ( U y) = n+(u)T(A)n+(uy (4) 

This equation merely expresses that T does not prefer any direction in Ti.. It would be a natural initial assumption 
for good cloning devices but, of course, in our case it will come out as a result of the minimization: T from equation 
is obviously covariant, because Sm commutes with all ir§ M (u). It is convenient to state the covariance condition 
as a fixed point property: we define the action r of unitary rotations on cloning maps by 

(t u T)(A) = ir+(uyT(n® A ^u)An® M (uyy+(u), (5) 

so that T is covariant iff t u (T) = T for all u £ U(o0. We denote by T the average of t u T with respect to u, i.e., 

T = J duT u (T), (6) 

where "du" denotes the normalized Haar measure on U(d). 

The fact that the cloning error A ono does not single out a direction on Ti either is expressed by the — easily verified 
— equation 

A one (r„T) = A onc (T). (7) 
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Similarly, we can get an estimate of A onc (T): The functional A onc is defined as the supremum of a set of convex 
expressions in T. Therefore, it is convex, and A ono (T) < A one (T). So as long as we are only interested in finding some 
cloning map with minimal A one , we may restrict attention to covariant ones. 

There is a similar simplification, which we can make "without loss of cloning quality" : A onc is invariant under 
a change of the ordering of the clones. That is to say, if V : Tl® M — ► 7^® M is a permutation operator, and if we 
define ryT by (t\/T)(A) = T(VAV*), we may replace T by its average over permutations without loss of cloning 
quality. That is, we may assume that ryT — T for all permutations V . We will refer to this property as permutation 
invariance. 

Our strategy is now to assume U(d)-covariance and permutation invariance of T, and to show that there is a unique 
solution to the variational problem with these additional properties. The above convexity argument then implies that 
no other cloning map can do better. But since the functi onal A one is not strictly convex, we will need an extra step to 
establish uniqueness. This we will do in subsection III F by showing that any cloning map whose mean is the optimal 
covariant doner has to be covariant itself. 



B. Reduction to the extremal covariant case 



The functional A onc involves only operators T(A) with A of the special form A = d( fe ) = 1} (g> a (g> 1®( M fc) e 

B(H® M ). Now due to permutation invariance T(a^)) does not depend on k, and we have 

T (««) = ^ T (E a «)- ( g ) 

k 

What makes this equation useful is that on the right hand side T is now applied to one of the generators of the 
representation n§ N : we have exp(i X)feLi a (k)) = (exp(ia))^ M . Because T is covariant, we can determine how the 
operators in equation @ transform under U(d)-rotations: 



N 



(u)T(a( k ))ir%(u)* = T((uau*) {k) ), (9) 



where the multiplication of a and u on the right hand side is in the d x <i-matrices. This property fixes the "transfor- 
mation behaviour" of the operators T(a^))> an d as we will see, this essentially fixes the tuple of operators T(a^))- 
Of course, a = 1 in (||) simply leads to T(l^)) = T(T) = 1. The operator il is the (anti-hermitian) generator of the 
subgroup of unitaries multiplying each vector with the same phase. More interesting are the generators of SU(d), in 
which such trivial phases have been eliminated. These generators, in other words the Lie algebra su(d), are exactly the 
traceless anti-hermitian d x d-matrices. In the qubit case (d — 2) su(2) is spanned by the Pauli matrices (multiplied by 
i), and 3-tuples of operators transforming like the generators are known in physics literature as "vector operators" . It 
is well-known that, due to the simple reducibility of SU(2), each irreducible representation of SU(2) contains exactly 
one vector operator (up to a factor), namely the generators (angular momentum operators) of the representation 
themselves. So all operators T(a^) are determined by the single numerical factor relating the operators T(a^) to 
the generators of the irreducible representation tt^. 

It turns out that the same idea works in the SU(ci)-case for arbitrary d. In order to state it precisely, we need a 
notation for the Lie algebra representation associated with a unitary representation of a Lie group. We define dn(X) 
to be the anti-hermitian generator of the one-parameter subgroup generated by X, i.e., 



Then the desired property of a representation is stated in the following definition: 



(10) 



2 Definition. Let ir : G — > B(Htt) be a finite dimensional unitary representation of a Lie group G with Lie algebra 
g. Then g is said to be non-degenerate in B(?t K ) with respect to tt, if any linear operator L : g — > B(l-t K ) with the 
covariance property ■K(g)L(X)n(g)* = L(gXg~ 1 ) is of the form L(X) = AdTr(X), for some factor A G C. 

As we argued above, su(2) is non-degenerate in every irreducible representation of SU(2). However, for d > 3 we 
can find representations containing degenerate copies of the generators, and we have to make sure that the special 
representations occurring in the present problem arc of the "good" kind. This is the content of the following Lemma, 
proved in Appendix A. 

3 Lemma. su(d) is non-degenerate in B(TL® N ) with respect to tt^. 
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4 Corollary. Let it : U(d) — > B{H^) be a unitary representation, and let T : B(H 7T ) — > B(H+ ) be a completely 

+ (u)T(A)k+( 



positive normalized and XJ (d) -covariant map, i.e. T(tt(u)Att(u)*) = tt~^ ! (u)T(A)tt~^ ! (u)* . Then there is a number ui(T) 



such that 

N 



T(dTr(a)) = u(T)J2*(k), 



k=l 



for every a £ B(H) with tr(a) = 0. 

Given u>(T) for 7r = tt^ m , we can compute the cloning error A one (T) as follows: given a £ £?(7i) with < a < 1, 
we can write a = al + a' with tra' = 0. Then 



1 M 
T(a (fe) ) = al+— t(J2<*{i) 



i=i 

N 



Z=l 



and with a' = a - al and a 



r( OW )=^(i-^)i+^(x> ( ; 



In any state ip £ TL we get 



where j{T) — jjLu(T) is the Black-Cow factor already mentioned in Section |l[ With 

sup \{ip,aip) - ^ ] = ^— — 

we get 

d- 1 



A one (T) = 



^ |i- 7(r)l • (ii) 



We remark that the largest possible lu(T), to be determined below, still makes the second term in the absolute value 
less than 1, so we could omit the absolute value signs. In any case, we will only seek to maximize u>(T) from now 
on, ignoring the possibility of w(T) > M/N, anticipating that it will be ruled out by the result of the maximization 
anyway. 

An important observation about the Corollary and formula ( |iT| ) is that ui is clearly an affine functional on the convex 
set of covariant cloning maps (i.e., lo respects convex combinations). Whereas we previously used the convexity of A ono 
to conclude that averaging over rotations and permutations (and hence a move towards the interior of the convex set 
of cloning maps) generally improves the cloning quality, we now see that the optimum can be sought, as for any affine 
functional, on the extreme boundary of the subset of covariant cloning maps. Therefore our next steps will be aimed 
at the determination of the extremal U(d)-covariant and permutation invariant cloning maps, and, subsequently the 
solution of the variational problem for these extremal cases. 



C. Convex decomposition of covariant cloning maps 

For the first reduction step we use the close connection between the permutation operators on 7-^® M and the 
representation 7r§ M . Let (tt® M ) denote the algebra of all operators on Ti® M commuting with all n§ M (u) = u® M . 
This algebra consists precisely of the linear combinations of permutation unitaries |l^, Theorem IX. 11. 5]. So consider 
a reduction of 7r|f M into irreducibles, i.e., an orthogonal decomposition of the identity into minimal projections E a £ 



5 



(ir® M )' . Then due to covariance the operators T(E a ) commute with all 7ri(u), and because the latter representation 
is irreducible, they must be multiples of the identity, T(E a ) = r a l, say. Because T(VA) = T(AV) for permutation 
operators V, we also have T(AE a ) — T(E a AE a ). Hence 

T a {A)=r- 1 T{E a AE a ) (12) 

is a legitimate cloning map in its own right (provided r a ^ 0). Moreover, 

T(A) — T(AE a ) = Y,T{E a AE a ) = ^r a T Q (A) (13) 

a a a 

is a convex decomposition of the given T into such summands. Maximizing lo(T) = ^ a r a u>(T a ) thus means con- 
centrating the coefficients r a on those a, for which oj{T a ) is maximal. At this stage it is perhaps already plausible 
that only the summand T a , for which E a = Sm is the projection onto the symmetric subspace will give the best 
u}(T a ), because this is the space supporting the pure states a® M the doner is supposed to approximate. In fact, for 
the optimization of A a n in Jl6| this idea leads directly to a simple solution. In the present case we found no direct 
proof of this plausible statement. 

We therefore have to enter into the further convex decomposition of each T a . The output states of this cloning 
map are supported by H a = E a TL® M , and we will restrict T a accordingly, i.e., we consider it as a covariant map 
T a : B(H a ) — > B(Hf N ) 7 which is covariant with respect to the restricted representation n a = 7r® M \H a and njj. 

As for any completely positive map, the convex decompositions of T a are governed by the Stinespring dilation [fL8| . 
Since we are looking, more specifically, for decompositions into covariant completely positive maps, we have to invoke 
a "covariant" version of the Stinespring dilation Theorem |]l^| , which is stated in Appendix B for the convenience 
of the reader. According to this Theorem we can write a covariant completely positive T a : B(H a ) — > B(Hf N ) as 
T a (A) = V*(A ® l/c)V, where K. is some auxiliary Hilbert space carrying a unitary representation tt : V(d) — > £>(/C), 
and V : 7i® N — > TL a ® K, is an isometry intertwining the respective representations, i.e., 

Vtt+(u) = (7r («)<8>7r(u))V. (14) 

The convex reduction theory of T a is now the same as the reduction theory of tt into irreducibles: if Fp is a minimal 
projection in the algebra w' , and hence ir\FplC is irreducible, then A ^ V* (A® Fp)V is a covariant map, which cannot 
be further decomposed into a sum of covariant completely positive maps (see Appendix B). Note that V*(l ® Fp)V 
commutes with the irreducible representation 7r^, so that once again this summand is normalized up to a factor: 
V*(1®F/3)V = r p l. Therefore T a = J^p rpT af3 , where each T af} (A) = r^V* (A® Fp)V is again an admissible cloning 
map. The following statement summarizes the result of the decomposition theory of T. 

5 Proposition. Let T : B(H.® M ) — » B(H® N ) be a V(d) -covariant and permutation invariant cloning map. Then T is 
a convex combination T — YlaB r apT a i3 such that each T a p is of the following special form: T a/ 3(A) = V*(A ® 1p)V, 
where V is an intertwining isometry between tt~^ and n a ® up, such that Tr a : V(d) — » B{TL a ) is an irreducible 
subrepresentation of ir§ M , and up : U(d) — > B{Hp) is also an irreducible unitary representation. 

This Proposition summarizes all that is needed for the further treatment of the variational problem. However, we 
could have made a slightly stronger statement by eliminating the non-uniqueness introduced by the choice of the 
minimal projections E a . If the subrepresentations 7r a and 7r a ' are unitarily equivalent, then they can be connected by 
a unitary, which is again a linear combination of permutations. Hence the contribution of the term r a T a — r a pT a p 
to ui(T) depends only on the isomorphism type of 7r a . 

What we cannot assert in general, however, is that V is determined by the isomorphism types of 7r a and Tip: among 
the groups SU(d) only d = 2 is "simply reducible", which means that the space of intertwiners between 7r 7 and 
7r Q ® 7Ty3 is at most one-dimensional for arbitrary irreducible representations n a , np , 7r 7 . In the following subsection 
we will therefore focus on the qubit case, and show how to determine oj(T a p) from the representations involved. This 
procedure will then be generalized to arbitrary d, and it will turn out that, perhaps surprisingly, in the general case 
^(Tap) also depends on ir a ,irp only up to unitary equivalence. 



D. Maximizing lo in the case d = 2 

For d = 2 the representations of SU(2) are conventionally labelled by their "total angular momentum" j — 
0,1/2,1,.... The irreducible representation 7Tj has dimension 2j + 1, and is isomorphic to with N — 2j in 
the notation used above. For j = 1 we get the 3-dimensional representation isomorphic to the rotation group, which is 



G 



responsible for the importance of this group in physics. In a suitable basis X\, X2, X3 of the Lie algebra su(2) we get 
the commutation relations [X%, X2] — X3, and cyclic permutations of the indices thereof. In the j = 1 representation 
diri(Xk) generates the rotations around the fc-axis in 3-space. The Casimir operator of SU(2) is the square of this vec- 
tor operator, i.e., C 2 = X)fc=i In the representation ttj it is the scalar j(j + 1), i.e., if we extend the representation 
dir of the Lie algebra to the universal enveloping algebra (which also contains polynomials in the generators), we get 
97Tj(C 2 ) = j(j + 1)1. We can use this to determine u}(T a p) for arbitrary irreducible representations. This computation 
can be seen as an elementary computation of a so-called 6j-symbol (see also for a context in which the same 
computation arises), but we will not need to invoke any of the 6j-machinery. 

So let V be an intertwining isometry between 7r 7 and ix a (g> 7173, where a,/3, 7 £ {0,1/2,...} label irreducible 
representations. Then to is defined by 

uj-d^(X k ) = V*(dTr a {X k )®l0)V. (15) 

We multiply this equation by dn 7 (X k ), use the intertwining property of V in the form Vd7r 7 (X) = (dTr a (X) 01^ + 
1 Q <8> dwj3(X))V, and sum over k to get 

w • <9tt 7 (C 2 ) = V*(dTT a (C 2 ) ® 1 P )V + V*(dw a {X k ) ® d7r f3 (X k ))V. (16) 

k 

The tensor product in the second summand can be re-expressed in terms of Casimir operators as 

^(57r Q (X fc ) ® dTr (X k )) = i ^(57r Q (X fe ) ® 1/3 + l a ® dn f3 (X k )) 2 - ia7r Q (C 2 ) ® 1^ - ® ^ Q (C 2 ). 

Inserting this into the previous equation, using the intertwining property once again, and inserting the appropriate 
scalars for <97r(C 2 ) = C 2 (7r)l, we find that ui ■ C 2 (7r 7 ) = C 2 (7r Q ) + ^(C' 2 (7r 7 ) — €2(^01) ~ ^2(^/3)), and hence 

. = i + ^i^M. (17) 

2 2C 2 (^ 7 ) 

Note that we have only used the fact that the Casimir operator C2 is some fixed quadratic expression in the generators. 
This is also true for SU(d). Hence equation ([l?]) also holds in the general case. In particular, we have shown that 
for the purpose of optimizing uj(T a p) only the isomorphism types of it a and 7Tg are relevant, but not the particular 
intertwiner V. 

Specializing again to the case d = 2, we find 

1 a(a + 1) - /3Q3 + 1) 

W= ; ; . 18 

2 27(7+1) k ' 

Here 7 = N/2 is fixed by the number N of input systems, a is constrained by the condition that 7r Q must be a 
subrepresentation of 7r|^f/ 2 , which is equivalent to a < M/2. Finally, /3 is constrained by the condition that there 
must be a non-zero intertwiner between 7r 7 and n a ® ftp. It is well-known that this condition is equivalent to the 
inequality \a — /3\ < 7 < a + (3. This is the same as the "triangle inequality": the sum of any two of a, /3, 7 is larger 
than the third. The area of admissible pairs (a, (3) is represented in Fig. 1. 



Since x i— > x(x + 1) is increasing for x > 0, we maximize ui with respect to (3 in equation (18) if we choose (3 as 



small as possible, i.e., /3 = \a — 7|. Then the numerator in equation (18) becomes 

a(a + 1) - (3({3 + 1) = 2^7 - -f 2 + max{7, 2a - 7}, 
which is strictly increasing in a. Hence the maximum 

M + 2 



(19) 



is attained for and only for a = M/2 and (3 — (M — N)/2. 
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Note that the seemingly simpler procedure of first maximizing a and then minimizing /3 to the smallest value 
consistent with a = M/2 leads to the same result, but is fallacious because it fails to rule out possibly larger values 
of lu in the lower triangle of the admissible region in Fig. 1. The same problem arises for higher d, and one has to be 
careful to find a maximization procedure which takes into account all constraints. 



E. Maximizing u in the general case 

Let us generalize now the previous discussion to arbitrary but finite d. In this case irreducible representations of 



U(d) are labelled, according to Section A. 2 by their highest weight m = (mi, . . . , ma)- Hence we can decompose 
T : B(H® M ) -> B(Hf N ) as described in Prop. 5 into the sum T = £ (m a]( - w r mn T m , n , taken over the set 

W = {(m, n) e 1\ x 7L% | 7r m C tt® m and tt+ C 7r m <g> 7r n }. 

Here is an abbreviation for the set of all possible weights of irreducible U(d) representations, i.e. = 
{(mi, . . . , m d ) mi > m 2 > ■ ■ ■ > m d }. 

Our task is now to determine (m,n) 6 W such that lu = u>(T mtTl ) becomes maximal. To this end we consider in 
analogy to ( |l5| ) the equation 

Lo-dir+(X) = V*(dir m (X)®l n )V, Vlesu(d) (20) 

where V is an intertwining isometry between 7r^ and ir m (g) 7r n . Note that equation (|2^) is valid only for X e Su(d) 
(and not for X 6 u(d) in general). Hence we have to consider the second order Casimir operator C2 of SU(d) which 
is given, according to Appendix [A 5| , by an expression of the form C2 = J2jk S^XjXk- This is all we needed in the 
derivation of equation ( |l7| ) in the SU(2)-case. The generalization to arbitrary d hence reads 

1 , C 2 (n m ) - C 2 (ir n ) 

LJ = - H ; . (21) 

2 2C 2 (^+) 



The concrete form of C2(7r m ) as a function of the weights m is given in equation (|Aq ), and will be needed only later. 
Since 62(71"^) is a positive constant we have to maximize the function 

W 3 (m, n) ^ F(m, n) = C 2 (^ m ) - C 2 (7r n ) € Z (22) 

on its domain W. 

The first step in this direction is to reexpress ^(m, n) in terms of the U(d) Casimir operators C 2 and C\ . Note in 
this context that although equation rt2G) is, as already stated, valid only for X G su(d) the representations 7r m and 



7r n are still U(<i) representations Hence we can apply the equation C 2 = C 2 — 4Cf given in Section A 5 
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F(m,n) = C 2 (TT m ) - C 2 (^ n ) - -{Cf{n m ) - C?(tt„)). 



(23) 



This rewriting is helpful, because the invariants Ci turn out to be independent of the variational parameters: Since 
TTm C ir® M , and 97r^ M (l,i) = Ml, we also have Ci (ir m ) = M. On the other hand, the existence of an intertwining 
isometry V with Vnjj = n m ® 7r n T^ implies 

V<7i(7r£)l = V^CCi) = {dir^id) ® I n + I m ® 07r n (Ci)) 7 = (d(7r m )l + d(^ n )l) V 

and therefore Ci(7r t f ) = Ci(7r m ) + Ci(7r n ). Since Ci(n^) — N and Ci(7r m ) = M we get Ci(7r n ) = N — M. Inserting 
this into equation fl23h we find the functional 



F(m, n) = i*i(m, n) 



2MAT - 



where only Fi depends on the variational parameters, and is expressed explicitly (see equation A7) as 



d d 

W 3 (m, n) i > fi(m, n) = C 2 (^ m ) - C 2 (^ n ) = ^)(m| - n|) + j^(d - 2fc + l)(m fe -n t )eZ 



(24) 



(25) 



fe=i 



which remains to be maximized over W . 

To do this we have to express the constraints defining the domain W more explicitly. We have already seen that 
m G U\_ has to satisfy the constraint Sj=i m j = ^ n addition we get, due to equation Al ma > 0. To fix the 
constraints for n note that according to equation ( A4 ) 7r^ C n m (£> 7r n is equivalent to ir m C ® ir^ . Here we have 
introduced n = (ni, . . • , n<j) = (— TI4, . . . , — ni) as a notation for the highest weight of the representation conjugate 
to 7r n (i.e. 7f^ = 7Tn). Now we can apply equation (A3) to get 



C 7r m ® 7r n n fc = rrik ~ Mfc witn < ^fc < m fe - m fe+ i Vfc = 1, . . . , d — 1 and ^ fc = AT. 

In other words 



fc=i 



= {(m, n) I n = m — /i, and (m, //) 6 Wi} 



with 



Wi = {(m, /iJeZfxZ^ra^M, ^ ^ k — N and < ^fe < Wfe - mfe+i Vfc = 1, . . . , d - 1}. 



fc=i 



fc=i 



The function Fi can now be re-expressed in terms the new variables (m, //). To this end note that C2(7r n ) = C2(7r n ) 
C2(7r„). Hence we have 

Ft (m, n) = Fx (m, S) = Fx (m, m - /x) 



and therefore with equation (|25|): 



F^m, n) = ^ Mfc(2TO fc - 2fc - fi k ) + (d + 1) ^ ^ fe = F 2 (m, fi) + (d+ 1)N 

k=l k=l 

with the new function 

d 

Wi 3 (m,fj,) i-> F 2 (m,/i) = ^^(2TO fe -2k- fik) € Z. 



fc=i 



Hence we have reduced our problem to the following Lemma: 



(26) 



(27) 







6 Lemma. The function F2 : W\ — *• Z defined in equation (|27[) attains its maximum for and only for 



= (M, 0, ... ,0) and jM 



max 



( (A, 0, . . . , 0) for A < M 

\ (M,0, ...,0,7V -M) forN>M. 



Proof: We consider a number of cases in each of which we apply a different strategy for increasing F 2 . In these 
procedures we consider d to be a variable parameter, too, because if fi d — m d = 0, the further optimization will be 
treated as a special case of the same problem with d reduced by one. 

Case A: fi d > 0, fii < rrii — mj+i for some i < d. 
In this case we apply the substitution fii i— > (fii + 1), fid (fJ-d — 1) 5 which leads to the change 

5F 2 = 2(— /if + Md + (d- i - 1) + m i+ x - m d ) > 2(fi d + (d- i - 1)) > 

in the target functional. In this way we proceed until either all fii with i < d satisfy the upper bound with equality 
(Case B below) or fi d = 0, i.e., Case C or Case D applies. 

Case B: fi d > 0, fii — m^ — TOj+i for all i < d. In this case all fik, including fi d are determined by the rrik and by 
the normalization {fid = A — mi + m d ). Inserting these values into F 2 , and using the normalization conditions, we 
get F 2 (m, n) = F 3 (m) - 2(M + diV) - TV 2 with 

F 3 (m) = 2(JV + d)mi 

constrained by mi > • • • > > , mfc = M, and mi — m^ < N. 

k 

This dehnes a variational problem in its own right. Any step increasing mi at the expense of some other rrik increases 
i*2- This process terminates either, when M — m\, and all other rrik — 0. This is surely the case for M < N, because 
then fid — A — mi + m d > N — M > 0. This is already the final result claimed in the Lemma. On the other hand, the 
process may terminate because fid reaches or would become negative. In the former case we get fid = 0, and hence 
Case C or Case D. The latter case (termination at fid = 1) may occur because the transformation mi 1— > (mi + 1), 
"m-d l— * (m<j — 1) changes fi d = N — mi + m d by —2. There are two basic situations in which changing both mi and 
m d is the only option for maximizing F3, namely d = 2 and mi = = ■ ■ ■ = m d . The first case is treated below 
as Case E. In the latter case we have 1 — N — mi + m^ = N . Then the overall variational problem in the Lemma is 
trivial, because only one term remains, and one only has to maximize the quantity 2m^ — 2k — 1, with trivial maximum 
at k = 1, mi = M. 

Case C: fid — 0, m d > 0. For fid = 0, the number m^ does not enter in the function F 2 . Therefore, the move m^ 
and mi 1— » mi + m d , increases F 2 by fiim d > 0. Note that this is always compatible with the constraints, and we end 
up in Case D. 

Case D: fid = 0, m d = 0, d > 2. Set d 1— > (d — 1). Note that we could now use the extra constraint fi d ' < rn d i, where 
d' = d— 1. We will not use it, so in principle we might get a larger maximum. However, since we do find a maximizer 
satisfying all constraints, we still get a valid maximum. 

Case E: d = 2, fix = mi — m 2 , fi% = 1. In this case m = (mi,m 2 ) is completely fixed by the constraints. We 
have: mi + m 2 = M and fix + fi 2 = mi — m 2 + 1 = N hence mi — m 2 = N — 1. This implies 2mi = M + N — 1, 
2m 2 = M - N + 1 and since m 2 > we get M > N - 1. If M = TV - 1 holds we get mi = iV - 1 = M, m 2 = and 
consequently fix = N — 1. Together with fi 2 — 1 = N — M these are exactly the parameters where F 2 should take its 
maximum according to the Lemma. Hence assume M > N. In this case fi 2 — 1 implies that F 2 becomes NM — 3iV — 4, 
which is, due to M > TV, strictly smaller than F 2 (M, 0; iV, 0) = 2AfiV - N 2 - 2N. 

Uniqueness: In all cases just discussed the manipulations described lead to a strict increase of F 2 (m, fi) as long 
as (m, fi) ^ (m max , ^ max ) holds. The only exception is Case C with fii = 0. In this situation there is a 1 < k < d 
with fik > 0. Hence we can apply the maps d 1— » d — 1 (Case D) and m^ 1— » and mi 1— » mi + m^ (Case C) until 
we get fi d 7^ (i.e. d reaches k). Since /^i = the corresponding (m, fi) is not equal to (m roax , Mmax). Therefore we 
can apply one of manipulations described in Case A, Case B or Case E which leads to a strict increase of F 2 (m,fi). 
This shows that F 2 (m,fi) < F 2 (m max , ^ max ) as long as (m, fi) ^ (m max , /x max ) holds. Consequently the maximum is 
unique. QED. 



With this result and the equations fl2l|), (|22|), (p4[), ( |26| ) and ( P7| ) we can easily calculate ui n 

~ M + d 

Umax = W i = 

and with (O) we get A(T) > A(f ) with A(f ) from Theorem 1. 
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F. Proving uniqueness 



One part of the uniqueness proof is already given above: there is only one optimal covariant cloning map, namely T. 
This follows easily from the uniqueness of the maximum found in Lemma 6 and from the fact that the representation 



njj is contained exactly once in the tensor product ttJ^ ® t^m-u ( see equation A3 and the discussion in Subsection 



ni q) . 

Suppose now that T is a non-covariant cloning map, which also attains the best value: A one (T) = A one (T). Then 
we may consider the average of T of T (see equation (^)), which is also optimal and, in addition, covariant. Therefore 
T = T. The uniqueness part of the proof thus follows immediately from the following proposition: 

7 Proposition. Each completely positive, unital map T : B(H® M ) — > B(H® ) satisfying the equation T = T equals 
f. 

Proof: We trace back this statement to the main theorem of . To this end note that T = T implies the equivalent 
equation for the preduals: 



where t u acts on T* by: 

Tu T*(a) = tt® m (u)*T* (n+(u)an+(uy) n§ M (u). 

Furthermore we know from the main theorem of jl6| that tr (CT® M T„((T® Ar )) < gKU is true for all pure states a € i3*(W) 
and that equality holds iff T = T. Consequently we have 



tr (a® M T u T{a® N ))) du = S - tr {a 9M %{a 9lf )) 
) d[M\ 



tT(a® M T*(a® N )\ =0 



\d[M] v " v ") d[M] v * v " d[M] 

Since the integral on the left hand site of this equation is taken over positive quantities the integrand has to vanish for 
all values of u G U(d). This implies tr {a® M T{a® N )) = |rpgj for all pure states a € B*(TL). However this is, according 

to fil, only possible if T = f. QED. 
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APPENDIX A: REPRESENTATIONS OF UNITARY GROUPS 



Throughout this paper many arguments from representation theory of unitary groups are used. In order to fix 
the notation and to state the most relevant theorems we will recall in this appendix some well known facts from 
representation theory of Lie groups. General references are the books of Barut and Raczka [|T), Zhelobenko ]22| and 
Simon @. 



1. The groups and their Lie algebras 



Let us consider first the group U(d) of all complex d x d unitary matrices. Its Lie algebra u(d) can be identified 
with the Lie algebra of all anti-hermitian d x d matrices. The exponential function is then given by the usual matrix 
exponential X i— > exp(iX). u(d) is a real Lie algebra. Hence we can consider its complexification u(d) ® C which 
coincides with the set of all d x d matrices and at the same time with the Lie algebra Ql(d, C) of the general linear 
group GL(d, C). In other words u(d) is a real form ofg[(d,C). A basis offlI(d, C) is given by the matrices Ejk = \j){k\. 
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The set of elements of U(d) with determinant one forms the subgroup SU(d) of U(d). Its Lie algebra su(d) is the 
subalgebra of u(d) consisting of the elements with zero trace. Hence the complexification su(d) <8> C of SXi(d) is the 
Lie algebra of trace-free matrices and coincides therefore with the Lie algebra sl(d, C) of the special linear group 
SL(d, C). As well as in the U(d) case this means that su(d) is a real from of sl(d,C). The matrices Ejk are no longer 
a basis for sl(d, C) since the Ejj are not trace free. Instead we have to consider Ejk, j ^ k and Hj = Ejj — Ej+ij+i, 
j = 1, . . . , d— 1. The difference between s[(d, C) and gl(d, C) is exactly the center of Ql(d, C), i.e. all complex multiples 
of the identity matrix. In other words we have gl(d, C) = sl(d, C) © CI. A similar result holds for the real forms: 
u(d) = su{d) © EI. 

The (real) span of all iEjj, j = 1, . . . , d is a subalgebra of u(d) which is maximal abelian, i.e. a Cartan subalgebra 
of u(d). We will denote it in the following by t(d) and its complexification by ic(d) C $l(d, C). The intersection of t(d) 
with su(d) results in a Cartan subalgebra st(d) of su(d). We will denote the complexification by stc(d). Again the two 
algebras 1(d) and si(d) differ by the center of u(d) i.e. 1(d) = si(d) © EI and tc(d) = stc(rf) CI in the complexified 
case. 



2. Representations 

Consider now a finite-dimensional^] representation tt : U(c?) — > GL(7V, C) of U(d). It is characterized uniquely by 
the corresponding representation dir : u(d) — > gt(iV, C) of its Lie algebra, i.e. we have 7r(exp(A)) = exp(<9-7r(X)). The 
representation 8tt can be extended by complex linearity to a representation of Ql(d, C) which we will denote by 8tt as 
well. Hence dir leads to a representation 7r of the group GL(d, C). Similar notations we will adopt for representations 
of SU(d) and SL(d,C). 

Assume now that 7r is an irreducible representation of GL(d, C). An infinitesimal weight of tt (or simply a weight 
in the following) is an element A of the dual of t£(c£) of ic(d) such that dir(X)x = X(X)x holds for all X £ ic(d) and 
for a nonvanishing x <E <C N . The linear subspace V\ C C N of all such x is called the weight subspace of the weight 
A. The set of weights of tt is not empty and, due to irreducibility, there is exactly one weight m, called the highest 
weight, such that dTv(Ejk)x = for all x in the weight subspace of m and for all j, k = 1, . . . , d with j < k. The 
representation tt is (up to unitary equivalence) uniquely determined by its highest weight. On the other hand the 
weight m is uniquely determined by its values m(Ejj) = rtij on the basis Ejj of ic(d). We will express this fact in 
the following as "m = (mi, ■ ■ ■ ,rrid) is the highest weight of the representation 7r". For each analytic representation 
of GL(c£, C) the mj are integers satisfying the inequalities mi > 11J2 > . . . > mj, and the converse is also true: each 
family of integers with this property defines the highest weight of an analytic, irreducible representation of GL(d, C). 

In a similar way we can define weights and highest weights for representations of the group SL(c?, C) as linear forms 
on the Cartan subalgebra sic(d). As in the GL(d, C)-case an irreducible representation tt of SL(d, C) is characterized 
uniquely by its highest weight m. However we can not evaluate m on the basis Ejj since these matrices are not trace 
free. One possibility is to consider an arbitrary extension of m to the algebra tc(<i) = stc(d) © CI. Obviously this 
extension is not unique. Therefore the values m(Ejj) — rrij are unique only up to an additive constant. To circumvent 
this problem we will use usually the normalization condition = 0. In this case the integer rrij corresponds to the 
number of boxes in the j th row of the Young tableau usually used to characterize the irreducible representation w. 
Another possibility to describe the weight m is to use the basis Hj of stc(d). We get a sequence of integers lj = m(Hj), 
j = 1, . . . , d — 1. They are related to the rrij by lj = mj — Wj+i- Each sequence li, . . . , ld-i defines the highest weight 
of an irreducible representation of SL(<i, C) iff the lj are positive integers. 

Finally consider the representation 7f conjugate to tt, i.e. Jf(u) = tt(u) . If tt is irreducible the same is true for ?F. Hence 
7r admits a highest weight which is given by (— md, —md-i, ■ ■ ■ , —mi). If tt is a SU(d) representation we can apply the 
normalization rrid = 0. Doing this as well for the conjugate representation we get (mi, mi — md-i, ■ ■ ■ , mi — m-2, 0). In 
terms of Young tableaus this corresponds to the usual rule to construct the tableau of the conjugate representation: 
Complete the Young tableau of tt to form a d x mi rectangle. The complementary tableau rotated by 180° is the 
Young tableau of W. 



2 AW representations in this paper are finite dimensional. 
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3. Tensor products of representations 



Consider now two finite dimensional irreducible representations 7r m , 7r n of U(d) with highest weights m, n. Their 
tensor product 7r m <g) 7r n is completely reducible. If r„ denotes the multiplicity of the irreducible representation 7r 
in 7r m ® 7r n then this means that 7r m ® 7r n = (J)^ r^Tr. Hence to decompose the representation 7r m ® 7r„ we have to 
compute the integer valued functions (m, n) i— > rv(m, n). There are several general schemes to do this (see e.g. |2^, 
Ch. XII]). However we are only interested in the foll owing special cases. The highest weight of the representation 
7Ti : ~U(d) 3 U U E GL(d, C) (denoted 7To in Section [II A ) is 1 = (1, 0, ... , 0). Consider the N— fold tensor product 
of this representation It can be decomposed as follows 



•! + ■■•+», 



r(mi, . . . ,md)ir r , 



(Al) 



denotes the irreducible representation with highest weight (mi, 



where 7T mil ... 

r(mi, . . . , rrid) are determined by the following recurrence relation: 

r(mi, . . . , m d ) = r(mi - 1, . . . , m d ) + r(mi,m 2 - 1, ... , m<j) + 



, m d ) 



The coefficients 



+ r(mi, .. .,md - 1). 



(A2) 



Consider now the A^— fold symmetric tensor product of 7Ti (denoted 7r^ in Section III A). It is irreducible with 



highest weight Nl = (N,0, . . . , 0) (hence 7r^ = 7Tjvi)- The tensor product of this representation with an arbitrary 
irreducible representation 7r m (with highest weight m = (mi, . . . , m d )) is 



TTNl <£> 7Tm = ^ 7T mi+Ml , 

»<(it + l^ m )o- ra Hl 



From this equat ion w e also get a condition for ttni to be contained in an arbitrary tensor product 7r r 
need in Section [HE: For arbitrary weights m, n, p we have 



C 7T„ <g> 7Tr 



7T„ G 7T D (g) 7T r 



(A3) 

1 7r n which we 
(A4) 



If two irreducible representations 7r m , 7r n of SU(d) are given we can characterize them, as described above, by their 
highest weights m = (mi, . . . , m d ) and n = (m, . . . , n d ) using the normalizations m d — and n d — 0. After applying 
the stated theorems to the tensor product of the corresponding U(d) representations we can restrict the summands 
in the resulting spectral decomposition back to SU(d), i.e. we renormalize the hcighest weigths (mi, . . . ,m d ) to the 
m d — case. 



4. Nondegeneracy of su(d) 



We are now ready to discuss the group theoretic part of the proof of our main theorem, i.e. Lemma 3 which we 
have only stated in Section III . According to Def. 2 we have to show that each linear operator A : su(d) 
the covariance property 



n% N with 



Tr+^A^KCg- 1 ) = AigXg- 1 ) 



(A5) 



is of the form A(X) 



Xdn^(X) with a constant factor A. Here nj^ is the irreducible representation of SU(g0 introduced 



in Section III A. (Hence we have 7r^ = ~kni using the notation introduced in Subsection A3 of this appendix.) 

To reformulate this statement note first that the map g ^ tt + 
representation of SU(cQ on the representation space <8> 7^2 



N {g) ■ 7r^(<7 1 ) can be interpreted as a unitary 
. In fact it is (unitarily equivalent to) the tensor 
g _1 E B(su(d)) is the adjoint representation of SU(d) this implies that each 
map X satisfying ( |A5| ) intertwines 7r^ (g) 7r^ and the adjoint representation Ad. Note second that the representation 
dn^ of the Lie algebra su(d) satisfies equation (A5) in an obvious way (with A = 1) hence we have to show that all 



product irjj <g> 7r^. Since SU(d) 9 g >—> g 



7r^ exactly once. 



such intertwiners are proportional, or in other words that Ad is contained in 7r^ 

Let us dis cuss now the tensor product tt^ (g> 
(see Section A2) and consequently the highest weight of its conjugate is (N, . . . , N, 0). We can apply now equation 



The irreducible representation 7r^ has highest weight (N, 0, . . . , 0) 



( |A3| ) which shows that the adjoint representation whose highest weight is (2, 1, ... , 1, 0) is contained in 7r^ ® ir 
exactly ones. This shows together with our previous discussion that su(d) is nondegenerate in Hf N with respect to 



' N 
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5. The Casimir invariants 



To each Lie algebra g we can associate its universal enveloping algebra ©. It is defined as the quotient of the full 
tensor algebra ©„ gNo Q® N with the two sided ideal 3 generated by X ® Y — Y <g> X — [X,Y], i.e. © is an associative 
algebra. The original Lie algebra g can be embedded in its envelopping algebra 0byg3XnI + 3e0. The Lie 
bracket is then simply given by [X, Y] = XY — YX. Moreover © is algebraically generated by g and 1. Hence each 
representation dir of g generates a unique representation dir of © simply by dir(Xi ■ ■ ■ Xf.) = dir(Xi) ■ ■ ■dTr(Xk). If 
dir is irreducible the induced representation dir is irreducible as well. 

We are interested not in the whole algebra but only in its center 3(©), i.e. the subalgebra consisting of all Z S © 
commuting with all elements of ©. The elements of 3(©) are called central elements or Casimir elements. If dir is 
a representation of © the representatives dir(Z) of Casimir elements commute with all other representatives dir(X). 
This implies for irreducible representations that all dir(Z) are multiples of the identity. 

Consider now the case g = gl(d, C). In this case we can identify the envelopping algebra © with the set of all left 
invariant differential operators on GL((i, C) (a similar statement is true for any Lie group). Of special interest for us 
are the Casimir elements belonging to operators of first and second order. Using the standard basis £!y of gl(d,C) 



introduced in Section A 1 they are given by 



Ci }2 I'-u and C 2 = E ikE kj . 

j=l j,k=l 

Of course C\ is as well of second order and it is linearly independent of C2. Hence each second order Casimir element 
of © is a linear combination of C2 and C\ . 

If dir is an irreducible representation of gl(c?, C) with highest weight (mi, . . . , mj) it induces, as described above, an 
irreducible representation dir of © and the images of c?7r(Ci) and dir(C 2 ) are multiples of the identity, i.e. <9tt(Ci) = 
Ci(n)l and dir(C 2 ) = C 2 (tt)1 with 

d d 

Ci{ir) — '^^'Tij and Ci(jr) = / ^m? + 2_,( m j — m k)- (^-6) 

j=l j=l j<k 

Let us discuss now the Casimir elements of SL(d, C). Since SL(d, C) is a subgroup of GL(d, C) its enveloping algebra 
& is a subalgebra of ©. However the corresponding Lie algebras differ only by the center of gl(d, C). Hence the center 
3(6) of 6 is a subalgebra of 3(©)- Since sl(d, C) is simple there is no first order Casimir element and there is only one 
second order Casimir element C2 which is therefore a linear combination C2 = C 2 + aC\ of and C 2 . Obviously 
the factor a is uniquely determined by the condition that the expression 

d I d 

C 2 {tt) = Ci(tt) + aC\ (tt) = + ^(rrij - m k ) + ot ^ m 3 | (A7) 

3 = 1 j<k V =1 

with dir{G2) = C?2(7r)l is invariant under the renormalization (mi, . . . , ma) 1— > (mi +//,..., + ji) with an arbitrary 
constant /1. Straightforward calculations show that a — — h- Hence we get C 2 = C 2 — \C\ and 



/ d d 

C 2 (ir) = - l(d - 1) - ^2 TOj-rrifc +d^2(rrij - m k ) J . (A8) 

\ J I .///• 3<k J 

Alternatively C2 can be expressed in terms of a basis (Xj)j of sl(d,C). In fact there is a symmetric second rank 
tensor g^Xj ® X k £ sl(d, C) ® sl(d, C) such that C2 coincides with the equivalence class of g jk in 6. In other words 
C 2 = 'Yl,j} z g^ k XjXk holds which leads to 

C^l^^g^diriX^diriXj) (A9) 

jk 

for an irreducible representation ir of SU(d). 
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APPENDIX B: STINESPRING THEOREM FOR COVARIANT CP-MAPS 



In this appendix we will state the covariant version of Stincspring's theorem |l9[ which we have used in the proof 
of Theorem 1. However, as in the rest of the paper, we will restrict the discussion to finite dimensional Hubert spaces 
(i.e. only cp-maps between finite von Neumann factors are considered). 

8 Theorem. Let G be a group with finite dimensional unitary representations 7Tj : G — > B(Hi) (i = 1,2), and 
T : B(J-L%) — > B(Hi) a completely positive map with the covariance property ir 1 {g)T{X)iT 1 {g)* = T(ir2(g)Xn2(g)*). 

1. Then there is another Unite dimensional unitary representation % : G — » B(Tt) and an intertwiner V : Hi — > 
H 2 ®H with Vwi(g) = ir 2 <g> ttV such that T(X) = V*(X 1)V holds. 

2. IfT = T a is a decomposition ofT in completely positive terms, there is a decomposition 1 = F a °^ ^ e 
identity operator on H into positive operators F a £ B{H) with [F a , n(g)} = such that T a (X) = V*(X®F a )V 

We only sketch the main ideas of the proof. The first step is Stinespring's theorem in its general form |L8| : There 
exists a representation r\ : B(H.2) — > B(IC) of the C*-algebra B(H.2) on a Hilbert space IC and a bounded operator 
V : Hi — > IC such that T(X) = V*r](X)V holds. Up to unitary equivalence there is exactly one such triple (IC, V, n) 
such that the vectors 7r(A)Vip G IC with ip S Hi and A G B{H-i) span K.. 

It is this uniqueness, from which the representation tt of G is constructed. Indeed, the objects V g = 
VTTi(g), and r) g (X) = rj(Tt2(g)XTt2(g)*) form a Stinespring dilation of the completely positive map T g (X) — 
iTi(g)*T('K2(g)X'!T2(g)*)'^i(g), which by covariance is equal to T. Hence by "uniqueness up to unitary equivalence" there 
is a unique unitary operator U g G B{K) such that V g = Vni(g) = U g V, and f] g {X) = i](ir2(g)XTT2(g)*) = U g n(X)U* . 
This can be simplified a bit further by the observation that according to the second equation the operators 
U g = r](TT2(g))*U g commute with all i](X). It is easy to see that the U g are a representation, and hence so is U: 
we have U g U h = rifa(g))*U g rifa(h)*)U h = V^2(g))*V^2(9)^(h)*n 2 (gy)U g U h = ^(g)* ^{g)^)* n 2 {gY)U gh = 
v(T2(gh)*)U gh = U gh . 

For a proof of part (1) we now only need to invoke the observation that all representations of B{H2) are of the form 
r\ ~ id ® 1 with K. = H2 ® H. (Here "~" denotes a unitary equivalence, which we will include as a factor in V. Since 
U g commutes with all i](X) = X ® I, it is of the form U g = 1 ® n(g), which proves the assertion. 

The second part of the theorem stated for a trivial group G = {e} is also known as the Radon-Nikodyn theorem 
coming with the Stinespring theorem. In general it asserts the existence of a partition of the identity operator on IC 
into operators F a commuting with all r)(X), giving the decomposition of T as T a = V*i](X)F a V. Again, we can 
write these as F a = 1 ® F a . Since the F a V are uniquely determined by the T a , it is easy to see that covariance of 
T a is equivalent to F a = n g F a n*. 
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